Interplay between Fe and Nd magnetism in NdFeAsO single crystals 
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The structural and magnetic phase transitions have been studied on NdFeAsO single crystals 
by neutron and x-ray diffraction complemented by resistivity and specific heat measurements. Two 
low-temperature phase transitions have been observed in addition to the tetragonal-to-orthorhombic 
transition at Ts ~ 142 K and the onset of antiferromagnetic (AFM) Fe order below Tm ~ 137 K. The 
Fe moments order AFM in the well-known stripe-like structure in the (ab) plane, but change from 
AFM to ferromagnetic (FM) arrangement along the c direction below T* ~ 15 K accompanied by 
the onset of Nd AFM order below TNd ~ 6 K with this same AFM configuration. The iron magnetic 
order-order transition in NdFeAsO accentuates the Nd-Fe interaction and the delicate balance of 
c-axis exchange couplings that results in AFM in LaFeAsO and FM in CeFeAsO and PrFeAsO. 

PACS numbers: 74.70.Xa, 75.25.-j, 61.50.Ks, 75.30.Kz 



Since the discovery of high-temperature superconduc- 
tivity in the fluorine-doped compound LaFeAsOi-ajF^; 
with a superconducting transition temperature of 
Tc = 26 K— , enormous efforts have been made to un- 
derstand the crystallographic, magnetic, and supercon- 
ducting properties of the i?FeAsO compounds {R = rare 
earth elements) and their doped superconducting deriva- 
tives with TcS up to ro55K2-28^ r^j^^ parent -RFeAsO 
compounds undergo a tetragonal-to-orthorhombic phase 
transition and show antiferromagnetic (AFM) order of Fe 
and R momentg^ii^"— . Although doping suppresses both 
the lattice distortion and Fe magnetic order, magnetism 
is believed to play a major role in the superconducting 
pairing mechanism in pnictide superconductors demon- 
strated, e. g., by the observation of a magnetic resonance 
in the superconducting statei^. Detailed knowledge of 
the magnetic order and crystallographic structure is key 
to gaining insight into the magnetic and magneto-elastic 
coupling, but has been hampered, so far, by the absence 
of sizeable i?FeAsO single crystals. 

Here we report on a neutron and x-ray diffraction 
study of recently available^S NdFeAsO single crystals 
complemented by resistivity and specific heat measure- 
ments. In prior neutron powder studies, the tetragonal- 
to-orthorhombic lattice distortion has been observed be- 
low Ts ^ 150 KiS followed by the onset of the well-known 
stripe-like AFM order of Fe below Tn -IdlK^^. The 
onset of the Nd AFM order was reported for temper- 
atures below TNd 2 K in neutron powder diffraction 
experiment o^'^d^'^^ . However, resistivity^! and muon- 
spin relaxatioriSS, studies indicate a low-temperature 
phase transition at 5-6 K. Our NdFeAsO single crys- 
tals show two low-temperature transitions in addition 
to the high-temperature transitions at Ts ^ 142 K and 
Tn ^ 137 K. Below T* ~ 15 K, the Fe moments undergo 
a strongly hysteretic magnetic transition from AFM ar- 
rangement to a ferromagnetic (FM) arrangement along 



the c direction but unchanged AFM stripe-like order in 
the (ab) plane. The Nd moments order AFM below 
TNd '^6K, coupled to the Fe magnetic order by a com- 
mon magnetic unit cell, indicating a complex interplay 
between Nd and Fe magnetism. 

NdFeAsO single crystals were grown out of NaAs flux 
as described in Ref.|20j. Batches with mm-sized single 
crystals and masses up to 7mg were carefully examined 
by wavelength dispersive spectroscopy (WDS) in a JEOL 
JXA-8200 Superprobe electron probe microanalyzer and 
x-ray powder diffraction. As in Ref.[2^ WDS measure- 
ments on different samples from the same batch show no 
noticeable compositional variation and confirmed, within 
the limits of the measurement, the atomic ratio of 1:1:1:1. 
No impurity phases were detected beyond the NaAs flux. 

Specific heat and resistivity were measured in a Quan- 
tum Design Physical Property Measurement System?^ 
by the relaxation method and the four-probe technique 
using Epotck H20E silver epoxy contacts on platinum 
leads. The homogeneity of the batches was examined by 
resistivity measurements of several single crystals show- 
ing similar behavior. The reproducibility and excellent 
homogeneity of the grown material were additionally 
confirmed by neutron diffraction measurements on poly- 
crystalline samples consisting of more than 300 crystals. 
They showed well-defined transition temperatures and 
Bragg peaks in excellent agreement with the single crys- 
tal measurements reported here. To check potential fer- 
romagnetism, the magnetization, M, was measured with 
increasing temperature from 5 to 20 K for eleven fields, 
H, between 0.1 and 5T and the differential susceptibil- 
ity was calculated at each temperature. Parallel to the 
c axis, the zero-field intercept of the M vs H line char- 
acterizing a FM component was essentially temperature 
independent giving an upper limit of 0.012 /is/formula 
unit. Perpendicular to the c axis the value was an order 
of magnitude smaller. 
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The lattice distortion was characterized by high- 
resolution x-ray diffraction measurements using a four- 
circle difFractometer and Cu-_fi'Q,i radiation from a rotat- 
ing anode x-ray source, selected by a Ge(l 1 1) monochro- 
mator. The plate-like single crystal was attached to a 
flat copper sample holder on the cold finger of a closed- 
cycle displex refrigerator. A mosaicity of ^ 0.03 deg was 
determined from full width at half maximum of the rock- 
ing curve of the (0 6) Bragg peak. The orthorhom- 
bic distortion, S, was determined as described in Ref.[23 
using ^ 0) scans through the tetragonal (1 1 7) Bragg 
peak which splits into the orthorhombic (2 7) and (0 2 7) 
Bragg peaks below Tg. 

The neutron diffraction experiments were carried out 
using the BT9 thermal neutron triple-axis spectrometer 
at NCNR, NIST. The crystal with a mass of 3mg was 
mounted on a thin aluminum post, oriented in the or- 
thorhombic (hOl) scattering plane, and placed in an ILL 
cryostat. Rocking scans were resolution limited, indi- 
cating the high quality of the sample. Neutrons with a 
fixed incident energy of 14.7 me V and a coUimation of 
40'-48'-sample-40'-open were used. Additional measure- 
ments were taken on a number of crystals heat treated 
in various ways to check the results, and the observed 
transitions were all consistent within uncertainties. De- 
scriptions are given in terms of the orthorhombic crystal- 
lographic (nuclear) unit cell present below T5. 

In Fig. [U we compare the four measurements as a 
function of temperature. The anomaly at T5 ~ 142 K 
in specific heat and resistivity measurements correlates 
with the phase transition from tetragonal P A/nmm to 
orthorhombic C mma. The orthorhombic lattice dis- 
tortion, 6, increases below T5 for decreasing tempera- 
ture without further anomalies. Figure[3l[a) shows the 
splitting of the tetragonal (2 2 0) Bragg peak into the or- 
thorhombic (4 0) and (0 4 0) Bragg peaks from our high- 
resolution neutron diffraction measurements and demon- 
strates the bulk nature of the lattice distortion and the 
homogeneity of the sample, through resolution-limited 
peak shapes and a uniform transition temperature. 

The second anomaly found in specific heat and resis- 
tivity data at T/v ~ 137 K corresponds to the onset of the 
well-known AFM stripe-like Fe orderiS. The intensity of 
the magnetic (1 3/2) Bragg peak increases smoothly 
below Ttv [see Fig.[ljd)]. The inset illustrates that this 
transition shows no hysteresis upon cooling and warming 
through Tjv consistent with a second order transition. 

At low temperature, a weak and broad anomaly in the 
resistivity can be associated with the disappearance of 
the magnetic (1 3/2) Bragg peak at T* ~ 15K signal- 
ing a change in the Fe magnetic order. This behavior is 
demonstrated in Fig.|3l(c) by rocking scans through the 
(1 3/2) peak at selected temperatures. At the same 
temperature, T*, new magnetic Bragg peaks occur as il- 
lustrated in Fig.[3i;b) by the (1 1) peak and in Fig.HJd) 
by the (10 0) peak. The intensity of these peaks increases 
with decreasing temperature as the Fe order changes, 
while below TNd 6 K the increase is much more pro- 
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FIG. 1: (Color online) (a) Temperature dependence of the 
specific heat, (b) Resistivity (open symbol, left axis) and 
its derivative (line, right axis), (c) Orthorhombic distortion, 
S, determined by x-ray diffraction as described in the text, 
(d) Neutron peak intensity for the magnetic Bragg peaks at 
(10 0) (left axis) and at (1 3/2) (right axis). The insets 
show expanded scales. 



nounced. At TNd, a- sharp peak is also observed in the 
specific heat [see inset of Fig.[IJa)] and the temperature 
derivative of the resistivity indicating the onset of the Nd 
magnetic order. 

Details of the temperature dependence for these low- 
temperature transitions are shown in Fig. [21 The tran- 
sition at T* 15K shows hysteretic behavior illustrated 
in Fig.[2l^b), through the temperature dependence of the 
magnetic (1 3/2) and (101) Bragg peaks. The strong 
hysteresis of approx. 3 K points to the first order charac- 
ter of this transition at T* . As illustrated in Fig. [51 the 
(1 0) and (1 1) Bragg peaks follow the same trend with 
(i) weak intensity between T* and Twd, and (ii) a strong 
but smooth and monotonic intensity increase below Twd 
without reaching saturation down to the lowest measured 
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FIG. 2: (Color online) Neutron peak intensity for the mag- 
netic Bragg peaks at (a) (10 0) and (b) at (101) (left axis) 
and (10 3/2) (right axis). The inset in (a) shows an expanded 
scale between TNd and T* . 
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FIG. 3: (Color online) (a) Neutron diffraction 9-29 scans 
through the tetragonal (2 2 0) nuclear Bragg peak at 
T — 158 K and its splitting into two orthorhombic Bragg 
peaks at T = 110 K measured with reduced coUimation of 10'- 
lO'-sample-lO'-lO'. Neutron diffraction rocking scans through 
the magnetic Bragg peaks at (b) (101) and (c) (1 3/2) at 
selected temperatures. 



temperature of 1.5 K. Both low-teraperature transitions 
are reminiscent of Nd2Cu04, where the Nii-Cu couphng 
leads to abrupt changes in the Cu spin structure along 
with a significant induced Nd moment, before the Nd 
spontaneously orders^^"— . It will be interesting to de- 
termine theoretically if a similar type of Fe-Nd coupling 
is the origin of this behavior in NdFeAsO. 

For the determination of the magnetic structures, a 
series of 5 crystallographic and 7, 5, and 15 magnetic 
Bragg peak intensities were measured by rocking scans 
at r =30K, 10 K, and 1.5 K, respectively, and corrected 
for instrumental resolution. The crystallographic Bragg 
peaks were fit using the published crystal structure^"^ to 
determine the scale factor as the only refinable parame- 
ter. The value of the Fe magnetic moment was then de- 
termined from the intensity of the magnetic Bragg peaks 
considering equal domain averaging and using this scale 
factor and the magnetic structure models given below. 

At T* < T = 30 K < Tat we find the same stripe-like 
AFM structure previously reportedi^. The Fe moments 
are coUinearly aligned along the a axis with AFM ar- 
rangement in the a and c directions and FM in the 
b direction. The magnetic unit cell for this magnetic 
structure is doubled along the c direction with respect to 
the orthorhombic crystallographic (nuclear) unit cell and 
yields the observed magnetic Bragg peaks at positions 
indexed by a propagation vector of (1 1/2). The de- 
termined value of 0.54(3) hb for the Fe moment is larger 
than the previously reported value of 0.25(7) hb^- The 
reason for this difference is unknown and points to differ- 
ences between polycrystalline and single crystal samples 
(potential small deviations in stoichiometry or strain). 



However, we note that the magnetic peaks associated 
with iron ordering were not reported in several other neu- 
tron powder scattering experiment o^'^i^^'^^ . 

The magnetic diffraction pattern at Twd < T = 10 K < 
T* can be indexed using a magnetic unit cell identical 
to the orthorhombic crystallographic (nuclear) unit cell. 
In particular, no doubling along the c direction is now 
required, signaling a change from AFM ordering in this 
direction above T* to FM ordering below T*. Together 
with the unchanged stripe-like order in the (ab) plane, 
this magnetic structure is simply described by a propaga- 
tion vector of (10 0). The intensity of all measured mag- 
netic Bragg peaks can still be fit by Fe moments alone 
with a reduced moment value of 0.32(3) /is collinearly 
aligned along the orthorhombic a axis. 

At T < TncI, the strongly enhanced intensity of the 
magnetic Bragg peaks indicates the onset of the Nd mag- 
netic order. The description of the magnetic diffrac- 
tion pattern requires both Fe and Nd moments to be 
ordered. The positions of the magnetic Bragg peaks are 
consistent with the reported low-temperature magnetic 
structurei^ii^. However, further measurements are nec- 
essary for a detailed analysis of the magnetic structure 
with its extensive number of parameters. The poten- 
tial overlapping of Bragg peaks from twin-domains with 
different orientations together with the small orthorhom- 
bicity require experiments with higher resolution or the 
usage of detwinned samples. 

Summarizing, below Tjv ~ 137 K, the Fe moments or- 
der AFM with the ubiquitous stripe-like arrangement 
in the orthorhombic (ab) plane. For T < T* 15 K 
the c axis configuration changes from AFM as ob- 
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served in LaFeAsO- to FM as found in CeFeAsO-^ and 
PrFeAsO^i^. This temperature dependent change in 
the Fe magnetic order in NdFeAsO is unique in the 
i?FeAsO series and demonstrates that the magnetic cou- 
phng in the c direction is at the border between AFM and 
FM. Further studies of -RFeAsO compounds with other R 
elements are necessary to determine whether this dehcate 
balance in the magnetic coupling is controlled by steric 
effects or the magneto-crystalline anisotropy of the R el- 
ement caused by crystal-electric field effects transferred 
through a complex coupling of the Nd and Fe magnetism. 
The proximity between the characteristic temperatures 



for the order-order transition of the Fe magnetic order at 
T* to the onset of the Nd magnetic order at T^^ ~6K 
together with a common magnetic unit cell at low tem- 
peratures suggests the second case as the more likely sce- 
nario. 
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